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It is shown that the basic observed properties of the gamma-ray bursts (GRBs) are accounted for 
if one assumes that the GRBs arise by blueshifting the emission radiation of hydrogen and helium 
generated during the last scattering epoch. The blueshift generator for a single GRB is a region with 
a nonconstant bang-time function ts(r) (described by a Lemaitre - Tolman (L-T) exact solution of 
Einstein’s equations) matched into a homogeneous and isotropic (Friedmann) background. Blueshift 
visible to the present observer arises only on those rays that are emitted radially in an L-T region. 
The paper presents three L-T models with different Big Bang profiles, adapted for the highest and 
the lowest end of the GRB frequency range. The models account for: (1) The observed frequency 
range of the GRBs; (2) Their limited duration; (3) The afterglows; (4) Their hypothetical collimation 
into narrow jets; (5) The large distances to their sources; (6) The multitude of the observed GRBs. 
Properties (2), (3) and (6) are accounted for only qualitatively. With a small correction of the 
parameters of the model, the implied perturbations of the CMB radiation will be consistent with 
those actually caused by the GRBs. A complete model of the Universe would consist of many L-T 
regions with different tsir) profiles, matched into the same Friedmann background. This paper is 
meant to be an initial exploration of the possibilities offered by models of this kind; the actual htting 
of all parameters to observational results requires fine-tuning of several interconnected variables and 
is left for a separate study. 
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I. MOTIVATION AND BACKGROUND 

In the Lemaitre [1] - Tolman [2] (L-T) models, in 
which the bang-time function tsir) is not everywhere 
constant, radial light rays emitted at the Big Bang (BB) 
can display infinite blueshifts {z = —1) to all later ob¬ 
servers. This happeirs when the radial rays are emitted 
at the generic points of the BB, at which dts/dr ^ 0 
[3-5]. On the other hand, gamma-ray bursts (GRBs) 
are observed and are believed to originate at large dis¬ 
tances from our Galaxy, up to several billion light years 
[6]. The question thus arises: could GRBs have been 
emitted during the last scatteriirg epoch, together with 
the relic radiation now observed as the cosmic microwave 
background (GMB), and then blueshifted to gamma-ray 
frequencies by the mechanism mentioned above? 

For the blueshift mechanism to work, the GRBs would 
have to originate in regions that emerged from a locally 
delayed BB.^ The relic radiation is emitted a finite time 
after the BB, at the last-scattering hypersurface (LSH), 
so the observed blueshift must be bounded from below, 
z > zls > where zls is the blueshift acquired be¬ 
tween the LSH and the present time. The main technical 
problem to solve is this: can zls be sufficiently near to —1 
that, with the free functions of the L-T model suitably 
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^ Regions where the BB occurred earlier than in the background 
generate shell crossing singularities [3, 7] in addition to blueshifts. 


chosen, the frequencies are blueshifted from the range of 
the emission spectra of hydrogen and helium (the only 
elements present in large amounts during last scattering) 
to the gamma-ray range observed today? The present 
paper answers this question in the positive - see Sec. IX. 

Section II provides the most basic information on the 
GRBs. Section III is an introduction to the L-T models, 
and Sec. IV provides information on light propagation in 
these models. Section V discusses the method of calcu¬ 
lating and the properties of redshift in the L-T models. 
It is shown there that nonradial rays emitted at the BB 
display infinite redshift to all observers, independent on 
whether dfs/dr at the BB is zero or not. In Sec. VI the 
definition of the extremum-redshift hypersurface (ERH) 
and the method of determining it are recalled. In Sec. 
VII it is shown how the Friedmann models follow as a 
limiting case of the L-T model, and geometric param¬ 
eters of the now-standard AGDM model are presented. 
Section VIII presents the background Friedmann model 
used in this paper, it has A = 0, fc = —0.4. 

Section IX presents three L-T models that account for 
the observed range of frequencies of the GRBs. Section 
X shows how the third model qualitatively explains the 
limited duration of the GRBs. A quantitative modelling 
of the duration would require a much higher numerical 
precision, but it is shown that the model contains a pa¬ 
rameter that can be adapted to the actually observed 
durations. Section XI shows how the same Model 3 ac¬ 
counts (qualitatively) for the afterglows of the GRBs. 

In Sec. XII, nonradial rays passing through the L-T 
region are discussed. It is shown that the angular diam- 
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eter of the L-T region as seen by the present observer 
in the sky would be ^ 2 degrees, and beyond that cone 
the CMB radiation would not be perturbed. Thus, with 
the parameters of the model suitably corrected, the im¬ 
plied perturbations will be hidden within the 1-degree 
cone of the present resolution of the GRB detectors. The 
implied pattern of the perturbations could be observa- 
tionally tested when the resolution improves. 

Section XIII explains how the models of Sec. IX can 
create the illusion of the collimation of the GRBs into 
narrow jets, even though they are emitted isotropically 
by the L-T regions. 

Section XIV shows how the models deal with the large 
distances to the sources of the GRBs. If the GRBs are 
really generated at last scattering, then the distances to 
them are even larger than currently believed. However, 
with inhomogeneities and blueshifts present, the redshift 
is not a monotonic function of distance, and thus fails 
to be a distance indicator, while the present estimates 
of distances to the GRB sources assume a homogeneous 
background model and use redshifts of the afterglows for 
calculating distances. As a by-product it emerges that 
local blueshifts arise also on nonradial rays, but they are 
not visible to the present observer: they are overcompen- 
sated by redshifts acquired earlier and later along the ray. 

Section XV shows that the models used in the pa¬ 
per would allow one to accommodate up to ^ 10^ GRB 
sources in the sky at the present time. Decreasing the 
diameter of the L-T region by a factor / (necessary any¬ 
way for other purposes) would increase this number by 

P- 

Improvements in the model needed to achieve a full 
quantitative fit are discussed in Sec. XVI. Conclusions 
are summarized in Sec. XVII. 


that all GRBs have afterglows, but some of them were not 
detected for technical reasons (like the afterglow lasting 
for too short a time) [14]. 

(4) They are probably focussed into narrow jets. 

(5) Nearly all GRBs come from very large distances, 
from over 10® to several billion light years. 

Currently, there is no universally accepted explanation 
of origins of the GRBs. There exist competing attempts 
at explanation by known astrophysical phenomena such 
as gravitational collapse to a black hole, a supernova ex¬ 
plosion or a collision of ultra-dense neutron stars [6]. 

The models presented in Sec. IX account quantita¬ 
tively for the frequency range in property (I), for prop¬ 
erty (5), qualitatively for properties (2) and (3), and are 
consistent with the hypothetical Property (4). Refer¬ 
ences to these properties will be marked by bullets 
Modelling individual GRBs with full quantitative agree¬ 
ment will require more precise (and much more time- 
consuming) numerical fitting of the ts(r) profiles to the 
time-profiles of the GRB frequencies. This is left for fu¬ 
ture research. The present paper is a proof of existence 
of L-T models that reproduce the basic GRB properties 
via suitably adjusted tsir) profiles. 


III. THE L-T MODELS 


The metric of the L-T models is: 


ds^ = dt^ - 


R , 


I -h 2E{r) 


dr^ - R {t, r) (dz?^ -|- sin dip ^), 

. 

where E{r) is an arbitrary function. The source in the 
Einstein equations is dust; its (geodesic) velocity field is 


II. BASIC FACTS ABOUT THE GRBS 

The amount of information on the GRBs is enormous 
[8]. In this preliminary study we do not try to interpret 
all the observational data that are available. Instead, we 
concentrate on the most characteristic properties of the 
GRBs to show how they follow from a simple combined 
Friedmann/L-T model. The following properties of the 
GRBs need to be accounted for [6]: 

(1) Their frequencies extend from Vjmin ~ 0.24 x 10^® 
Hz to ^^niax ~ 1-25 X 10^® Hz [9]^ (see also Ref. [12]). 

(2) They typically last from less than a second to a 
few minutes [6], but a few examples are known of GRBs 
lasting from over two to about 30 hours [13]. 

(3) For most GRBs, longer-lived and fainter “after¬ 
glows” at larger wavelengths are observed. It is believed 


= (5“o. 


(3.2) 


^ Converted from keV to Hz hy u = E/h, where 1 keV = 1.6 X 
10“^® J [10] and h = 6.626 X 10“®^ J s [11]. The lowest energy 
was read out from Figs. 1 and 2 in Ref. [9] as 10 keV. 


Because of the property p = 0, this model is inadequate 
for describing times before the LSH. 

The function R{t, r) is determined by 

rP = 2E{r)+2M{r)/R+]-KR^, (3.3) 

M(r) being another arbitrary function and A being the 
cosmological constant. We consider models with R,t > 0, 
E > 0 and A = 0. The solution of (3.3) is then: 

M 

= —(cosh 7?-I), 

sinh?7-7; = — — — [t-tB{r)], (3.4) 

where Ib (r) is one more arbitrary function; the BB occurs 
a.t t = tB{r). The mass density is 


up = 


2M^r 

R^Rr' 


(^f SttG 

“ n2 ■ 


(3.5) 
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The r-coordinate is chosen so that [5] 

M = Mor^, (3.6) 

and Mq = 1 (kept in formulae for dimensional clarity). 

The units used in numerical calculations were intro¬ 
duced and justified in Ref. [15]. Taking [10] 

1 pc = 3.086 X 10^^ km, 1 y = 3.156 x 10^ s, (3.7) 

the numerical length unit (NLU) and the numerical time 
unit (NTU) are defined as follows: 

1 NTU = 1 NLU = 9.8 x 10^° y 

= 3 X 10^ Mpc. (3.8) 


The L-T models are generalisations of (not alterna¬ 
tives to!) the Friedmann models used in astrophysical 
cosmology to describe our actual Universe. The relation 
between these two classes is explained in Sec. VII. 

The following equations that hold in an L-T model 
with E 0 [7] will be useful further on: 


R 


ir 


R,tr 


( M,r 

V M 


E,r 

E 


R 


(3.9) 


3 E,r M,r ^ , 


M 


R 




M 


i?,. 


(3.10) 


3 E,r 
2^ 


M 


{t — ts) — tB,r 


IV. LIGHT RAYS IN AN L-T MODEL 


The general equations defining the tangent vectors 
k°‘ = da;“/dA to geodesics of the metric (3.1), with A 
being the affine parameter, are 


dl + TW > 


-l- RR,t 


[k^)^ + sm^d {k'^y 


= 0 , 


(4.1) 


^ + 2^ k^k^ 
dX R,r 

(1 -H 2E)R 


R, 


E,., 


R, 


R,r 1 + 2E 

2 . . 9 „ ,on2 


{k^y 


(A’’) -hsin^i? {k^f’f =0, (4.2) 

Tt,r i- -I 

^3— -I- 2^^ k*k'^ + 2^^ k^k"^ — cos-dsiii-d {k'^y = 0, 
dA R R 

(4.3) 

^ + 2 ^ A‘A- + 2^ A^A- + 2^ A’^A- = 0. 
dA it it sinv 

(4.4) 

The geodesics determined by (4.1) - (4.4) are null when 


(A‘) - 


2 R,ym[_j^2 


1 + 2E 


(A’^) -|-sin^'d (A”^)^ =0. (4.5) 


Using R,t A* -I- R,r k'' = di?/dA and recalling that k^ = 
di?/dA, the general solution of (4.4) is 

R^sin^dA^ = Jo, (4.6) 


where Jq is constant along the geodesic. The special case 
Jq = 0 corresponds to two situations: 

(a) k"^ = 0, i. e. being constant along the ray, with 
•d being, as yet, unspecified, or 

(b) i9 = 0, i.e. the ray proceeding along the axis of 
symmetry, with (f being undetermined. 

Using the above, the general solution of (4.3) is 

R^ (k^y sin^ d + Jo^ = C^sm^i9, (4.7) 

where is another constant along the geodesic. When 
C = 0, the geodesic is radial. Then Jq = 0 and either 
(a) ■(? = 0 with (p being undetermined or (b) i9 has any 
constant value along the ray and (p is constant in conse¬ 
quence of (4.6). When C = ±Jo ^ 0, the geodesic re¬ 
mains in the equatorial plane i? = -k 12. Along any single 
geodesic, 1 ? = 7r/2 can be achieved by a transformation 
of the {t),(p) coordinates. 

For rays with Jq y 0, eq. (4.6) implies in addition: 

k'^ = ^ —>• 00 when i? —0. (4.8) 

dA 

Thus, if such a ray has |dr/dA| < 00 at the intersection 
with the BB, then dip/dr —>■ 00 , i.e. each of these rays 

meets the BB being tangent to a surface of constant r. 
From (4.6) and (4.7) we get 


{k^y + ^ {k^y = 


and then (4.5) becomes 


{^y 


R,y{Ry 

1 + 2E • 


(4.9) 

(4.10) 


Using (4.6) - (4.10), eqs. (4.1) - (4.4) simplify to 


dt 

dA 

dA* 

dA 





R,r 


C^R,t 

i?3 


dr 


^ Vl + 2E 
R,r 



C2 


AU 


(4.11) 

(4.12) 

(4.13) 

(4.14) 


The initial data for (4.11) - (4.14) are (t,r) = {to,To) 
- the coordinates of the observation point. The values 
of A will not appear in the numerical calculations, in all 
graphs the independent variable will be r, so no value for 
X{to) has to be assumed. 

The sign in (4.13) is -I- on those segments of the rays, 
on which r is increasing, and — where r is decreasing. 
Care must be taken in numerical calculations at those 
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points where r changes from increasing to decreasing or 
vice versa; there fc’' goes through zero and changes sign. 

Since the affine parameter is defined up to a constant 
factor, one more initial condition is usually assumed; it 
can be achieved by rescaling A: 

k\to) = ±1 (4.15) 

(+ for future-directed, — for past-directed rays). Why 
this is convenient will be seen in Sec. V. 

With (4.15) we have from (4.10), 

<RHto,ro) = (4.16) 

the equality occurs when ko^ = 0, i.e. when the ray is 
tangent to an r = constant sphere at (t,r) = {to,ro). 

Given an initial point {to, To) and an initial direction 
coded in C, (4.11) - (4.14) determine fc* and k^ all along 
the ray.^ To find k"^ and k’^, one has to specify Jq, then 
solve (4.7) to find i?(A), and finally calculate k''^{X) from 

(4.6) . Note that, given C, there is a whole bundle of rays 
(labelled by Jq) that have the same A:*(A) and kX{\). 

In the graphs, only the rays lying in the subspace = 
7r/2 (on which Jq = iC) will be shown. To calculate 
them, in addition to (4.11) - (4.14) one has to integrate 

(4.6) , which becomes 


dA 


(4.17) 


with the initial value = 0. The “—rays” are mirror- 
images of those with -I-, and will not appear in the graphs. 

On past-directed radial rays, on which C = 0, using 
(4.10), the equations to be integrated numerically are: 


dt 

dA 

= k\ 

(4.18) 

dfc‘ 

dA 

II 

1 

to 

^ 7 

(4.19) 

fc” 

_ V1 + 2F; j 

it^r 

(4.20) 

dr 

dA 

= k\ 

(4.21) 


with the initial condition (4.15). The sign in (4.20) is 
e = -|-1 on past-inward and future-outward rays and — 1 
on other rays. 


® The initial values of and ip are irrelevant for (4.11) — (4.14) 
in consequence of (4.9): all rays originating on the same sphere 
(t, r) = {to, To) with the same value of C have fc* and fc’’ expressed 
by the same formulae. The value of is needed if we want to 
know k'^ and k‘^ as well. Even then, po is not needed because 
of spherical symmetry. It becomes needed when the equations 
A:“ =da:“/dA are integrated to find the path of the ray. 


V. THE REDSHIFT 


The general formula for redshift is [16] 


1 I - {Upk^^)e 
{uuk'^)o ’ 


(5.1) 


where k^ is an affinely parametrised vector field tangent 
to a ray connecting the light source and the observer, 
both comoving with the cosmic medium. The subscript 
“e” means “at the emission event”, “o” means “at the 
observation event”. 

Consider a ray proceeding from a spacetime point Pi to 
P 2 and then from P 2 to P3. Denote the redshifts acquired 
in the intervals [^ 1 ,^ 2 ], [^ 2 ,^ 3 ] and [PijPs] by 2 : 12 , 223 
and Zi 3 , respectively. Then 


1 + Zi3 — (1 -|- Z 12 ) (1 -|- 2 : 23 ) ■ (5-2) 


On past-directed rays, with being given by (3.2) 
and using (4.15), we have 


1 -|- 2 : = —fc*. (5-3) 

On future-directed rays, (4.15) and (5.3) change to 

k*{to) = +1) 1 -|- 2 ; = fc*. (5-4) 

The conditions (5.3) and (5.4) are compatible with (5.2). 

For nonradial rays, on which C yf 0, the last term in 
(4.10) will go to infinity when i? —>■ 0. Thus, at the BB, 
independently of whether the second term in (4.10) stays 
finite or becomes infinite 

lim fc* = lim 2 ; = 00 . (5.5) 

R —^0 R —^0 

Now imagine a bundle of rays originating at an off- 

center event Pq (to,ro,'i9o, (po), and going off from Pq 
to the past in all directions. Within this bundle, there 
will be a ray going radially away from the center and 
another one going radially toward the center; for both of 
them (7 = 0. Calculating the redshift in a vicinity of the 
BB for these two rays is a complicated thing, but power 
expansions of the redshift formula in P [4] and numerical 
integrations [5] both indicate that, for radial rays, 2 : —>• 
— 1 as P —t 0 when dt^/dr yf 0 at the intersection of the 
ray with the BB, and z —t 00 when dts/dr = 0. The 
value z = —1 implies that the observed frequency would 
be infinite (which means: very large for rays emitted close 
to, but later than the BB). The property 2 < 0 goes by 
the name blueshift, and 2 = — 1 means that the blueshift 
is infinite. Note that by the very definition of redshift, 

l-b2>0 (5.6) 


must always hold. 
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VI. THE EXTREMUM-REDSHIFT 
HYPERSURFACE 

Along a radial ray, dr/dA = fc’’ 0 (except possibly 
at the BB, where k* —>■ 0, see (4.20)). Therefore, using 
(5.3) and (4.20), Eq. (4.19) may be written as [17], [7] 

1 dz R-itr fa 

- - — = -e , =. ( 6 . 1 ) 

1 + 2 dr y/l + 2E 

Thus, R,tr = 0 is the locus of extrema of z along radial 
rays, called the extremum-redshift hypersurface (ERH). 

Further in this paper, we will consider an L-T model 
with 2E = —kr^, where k = constant. The way of deter¬ 
mining the ERH in such a model was described in Ref. 
[5], and we only copy the results. The value of ry (defined 
in (3.4)) on the ERH is determined by 


where 


X sinh^(r7/2). (6.3) 

Having found (numerically) a;(r), and thus also ri{r) from 
(6.3), we find t{r) on the ERH from (3.4): 

iERH(r) = tB{r) + {sinh[ 77 (r)] - 77 (r)}. (6.4) 

The ERH does not exist along those rays that hit the 
BB where = 0 and is not determined at r = 0 [5], but 
the limits at r —>■ 0 and at tB,r —>■ 0 of the solution found 
at r X tB,r ^ 0 may exist. In particular, (6.2) - (6.4) 
imply that ry —>■ 0 (i.e. fERH —>■ ts) at all points where 
tB,r 0, and also at r —> 0 if limr_>.o tB,r is finite. This 
means that on approaching all those points, the ERH and 
the BB become arbitrarily close to each other. 

Equation (6.2) makes no reference to the initial point 
of the geodesic arc. Consequently, the ERH is observer- 
independent. The extremum value of redshift will depend 
on the initial point, but the location of the extremum will 
not: the extremum of z along a given geodesic will occur 
always at the same r. 

Consider a radial ray proceeding to the past from an 
initial point that lies later than the ERH. The redshift 
on it increases from 0 to a local maximum, achieved at 
the ERH. Further down the ray, z initially decreases.^ If 
the ray could continue to the BB, 2 ; would either decrease 
to —1 (if dts/dr 7 ^ 0 at the intersection) or increase to 
infinity (if dts/dr = 0 there). However, the L-T model 
does not apply at times earlier than the LSH. Can z 


^ In the models considered further on, the ERH will have the topol¬ 
ogy of a thick-walled tea cup, and the radial rays will intersect 
it up to four times. The 2 along them will thus have up to two 
local maxima and up to two local minima. 


become, before the ray crosses the LSH, sufficiently neg¬ 
ative to shift the optical frequencies to the gamma-ray 
range? It is shown in Sec. IX that this is possible when 
the functions E{r) and tB{r) are suitably chosen, and the 
observer is put in the right spacetime region. 


VH. THE FRIEDMANN LIMIT OF THE L-T 
MODEL, THE ACDM MODEL AND THE 
LAST-SCATTERING INSTANT 


The Friedmann limit of (3.1) follows when and 
E/M'^/^ are constant. With the coordinate choice (3.6) 
this means 2E = —kr^, where k is the Friedmann curva¬ 
ture index. Then (3.4) imply R = rS{t), and 


ds'^ = dt'^ - S^{t) 


1 


I — kr'^ 


dr ^ + (d'd^ + sin'^ 1 ? dip^) 


while (3.4) reduce to 

S{t) = 
sinh rj — Tj = 


(7.1) 


Mo 


i-k) 
(-fc)3/2 
Mo 


(coshry — 1 ), 
(t - Ib) ■ 


(7.2) 


Since the Friedmann models are isotropic around every 
observer world line, every null geodesic is radial, so ( 6 . 1 ) 
can be used. It is then easily integrated to give 


l + z = S{to)/S{E), (7.3) 


where to and t^ are the instants of, respectively, the ob¬ 
servation and emission of the light ray (note that (7.3) 
trivially obeys (5.2)). 

The ACDM model, now most often used as the “stan¬ 
dard” cosmological model, is a solution of Einstein’s 
equations for the metric (7.1) with dust source and 
A: = 0 > A [15]: 



2 


[t — Iba) , 


(7.4) 


where t = tBA is the instant of the BB. It is characterised 
by the Hubble parameter. 


Ho = S,t /S\,^. 


t^to 


(7.5) 


and two dimensionless constants: the density parameter 
and the cosmological constant parameter 




1 / SttGpo 

iHo^ V ’ 


-A 


(7.6) 


that obey Q.m + AIa = 1 ; Po is the present mean mass 
density in the Universe. The Hubble parameter Elo in 
(7.5) is related to the Hubble constant Ro by 


Ho = Hole. 


(7.7) 
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The current observations imply [18] 

(a™,riA) = (0.32,0.68), -Ho = 67.1km/(s x Mpc). 

(7.8) 

The above, via (3.7) - (3.8), (7.6) and (7.7), leads to 

[i7o,-A] = [6.71 (NLU)-\91.849164 (NLU)“^] . (7.9) 

The age of the ACDM Universe is [18] 

r= 13.819 X 10® y = 0.141 NTU. (7.10) 

The end of the recombination epoch (the last scatter¬ 
ing instant) occurs when the temperature of the cosmic 
matter drops below the one needed for ionising the hydro¬ 
gen atoms. It will be assumed that this temperature is 
uniquely determined by the local mass density, pLS- Con¬ 
sequently, it will be the same along every matter world 
line and in every L-T model. At larger densities, the 
L-T and ACDM models do not apply. 

In the ACDM model the last scattering occurs at the 
redshift [18, 19] 

ZLS = 1090. (7.11) 

Using this value and (7.4), one can calculate the corre¬ 
sponding instant tns in the ACDM model from [7] 

l + 2LS = ^(T)/5(tLs), (7.12) 

where T is given by (7.10). Namely, using (7.9) and (7.10) 
in (7.4) one finds 

S{T) = 0.51743812113024401. (7.13) 

Next, using (7.4) and (7.12) one finds 

tLS = In (y + + l) , (7.14) 

where 

Y (1 -h ZLs)“^/^ Sinh . (7.15) 

Then, using (7.9), (7.10) and (7.11), one finds® 

tLS = 4.86905016470083480 x 10"® NTU 

= 4.7716691614068183 X 10® y. (7.16) 

The value of S{t) att = t^s is now calculated from (7.4): 

S'ls = 0.000474278754473001721. (7.17) 

From (7.9), (7.8) and (7.6) using c = 3 x 10® m/s and 
(3.7) - (3.8) we obtain for the present mean mass density® 

npo = Sn^ino/cf = 43.223136 (NLU)”®. (7.18) 

Then, the density at last scattering is 

kpls = Kpo[S{T)/SLsf (7.19) 

= 56.1294161975316 X 10® (NLU)“®. (7.20) 


® See Appendix A for a clarification of the possible confusion con¬ 
nected with zls vs. the age of the Universe at last scattering. 

® This value is corrected with respect to Ref. [20]. 


VIII. THE BACKGROUND MODEL 

Each cosmological model in this paper will consist of 
a Friedmann background with an L-T island matched 
into it. The background model will be different from 
ACDM. For this preliminary study it is preferable to set 
the cosmological constant zero in order to be able to do 
exact calculations as much as possible. Our Friedmann 


background will have the following parameters: 


A 

= 0, 

(8.1) 

k 

= -0.4, 

(8.2) 

ts 

tBf = -0.13945554689046649 NTU 



~ —13.67 X 10® years. 

(8.3) 


The tef is taken from Ref. [5], where it was the asymp¬ 
totic value of the function tsir) in an L-T model that 
mimicked accelerating expansion (presented in Ref. [15]); 
it differs from (—T) given in (7.10) by ~ 1.6%. The value 
of k emerged in numerical experiments. 

For comparing the predictions of our models with the 
CMB data, we will need the redshift at last scattering 
in this background model. The density at last scattering 
Pls given by (7.20) must be the same as before, so S'ls 
must be the same, too. This time, however, the present 
value of S must be calculated for the background model 
with the parameters (8.1) - (8.3). It is, using (7.2) 

Snow = 0.45180345033414671. (8.4) 

The resulting redshift between the LSH and t = 0 in the 
background model, Snow/SpS) comes out to be 

1 + zbg = 952.611615159. (8.5) 

This differs from (7.11) by ~ 12.7%. The present tem¬ 
perature of the CMB radiation is directly measured, so if 
(8.5) were taken for real, it would imply that the temper¬ 
ature of the background radiation at emission was higher 
by 12.7% than current knowledge tells us (i.e. ^ 3380 K 
instead of ~ 3000 K). But a more appealing way to cure 
this discrepancy would be to change ts or fe (or both) 
so as to make Snow larger. This would require increasing 
Its I or increasing \k\ (to produce larger S at a given t—ts 
with k < 0, larger |fc| is needed, see Fig. 17.1 in Ref. [7]). 
However, these changes would have to be accompanied 
by changes in other parameters of the models presented 
further on, and such changes would require laborious re¬ 
calculations. Therefore, for the beginning, we will show 
how the present model deals with some properties of the 
GRBs exactly, and with others only qualitatively, just 
to explore the existing possibilities. See Sec. XVI for a 
discussion of improvements. 

IX. FITTING THE L-T MODEL TO THE GRB 

FREQUENCIES (PROPERTY (1) IN SEC. H) 

As already explained, we assume that the last scatter¬ 
ing in any model occurs at the density given by (7.20). 
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The density along a ray is calculated using (3.5), and the 
value of 2 ; at the moment when p = pLs emerges from 
(5.3) during numerical integration of (4.18) - (4.21). 

The spectra of the GRBs do not have the black- 
body forms [21], so, if the GRBs arise by the mecha¬ 
nism described below, then different frequencies must be 
blueshifted independently; see Appendix B. 

As the cosmological model we take a Friedmann back¬ 
ground, with k and Ib given by (8.2) and (8.3), into which 
several L-T regions are matched. Each L-T region has a 
nonconstant tsir), and, being spherically symmetric, de¬ 
fines its own radial directions, independent of the other 
L-T regions. The BB in such a model consists of flat 
parts that give rise to the Friedmann background, and 
spherical humps of delayed BB that evolve into the L-T 
regions. In this section we consider three examples of sin¬ 
gle L-T regions surrounded by a Friedmann spacetime. 
The presence of other L-T regions would not perturb 
the evolution of any single one, but of course would in¬ 
fluence the propagation of light passing through them. 
The calculations presented here assume that there is no 
such intervening L-T region between the one that emit¬ 
ted the ray and the observer sitting in the Friedmann 
background. The perturbations caused by intervening 
L-T regions may be investigated separately. The central 
theme of this paper is the existence of the mechanism 
producing sufficiently strong blueshifts. 

We assume that the gamma rays in the GRBs originate 
as emission radiation of hydrogen and helium, the only 
elements present in large amounts during the recombi¬ 
nation epoch, together with the radiation that will later 
become the GMB. The rays emitted in the Friedmann 
background and rays emitted non-radially in the L-T re¬ 
gion become redshifted and evolve into the CMB. The 
rays that are emitted radially in the L-T region become 
Wweshifted and become the GRBs. 

The emission frequencies of hydrogen lie between 

= 4.054 X 10^3 Hz, (9.1) 

corresponding to the wavelength of 7400 nm, and 

i^Hmax = 3.2 X lO^"^ Hz, (9.2) 

corresponding to 93.782 nm [22], with the frequency of 
the most intense line being 

= 2.1876 X 10^^ Hz, (9.3) 

corresponding to 656.2852 nm. The most intense helium 
emission lines have wavelengths between 388 nm and 846 
nm, i.e. within the range of the hydrogen spectrum [23]. 

The maximum intensity of the CMB radiation by to¬ 
day is at the frequency of ^ 200 GHz [24]. Using (7.11) 
and Wien’s law it is easy to calculate that the frequency 
of maximum intensity at the time of emission must have 
been 2.18 x lO^'^ Hz. This falls between t'Hmin and J^Hmax- 

The z needed to shift the J^Hmin to the lowest observed 
frequency of the gamma-ray bursts [9], 

~ 0.24 X 10^^ Hz, (9.4) 


is 


l + -2max~ 1.689 X 10-^ (9.5) 

The z needed to shift the t'Hmax to the maximum 
recorded cosmic gamma-ray frequency [9], 

« 1.25 X 10^3 Hz, (9.6) 

is 

1 + .Zmin ~ 2.56 X 10“®. (9.7) 

Consequently, 2 : in the GRB models has to obey 


^min Z <C. Zjjiax' 


(9.8) 


Actually, if a model that accounts for a given 1 -|- z is al¬ 
ready constructed, then there is no problem with making 
1 -I- z larger; the challenge is to account for a sufficiently 
small l-|-z. Therefore, we will aim at constructing models 
that have z < Zmin and z < Zmax at the LSH. 

The notation may be confusing here: the Zmax (so de¬ 
noted because it is the greatest value of z that will be 
needed in GRB models) is associated with the minimum 
values of the emission- and GRB frequencies; similar con¬ 
fusion exists for the Zmin- 

The shape of the hump in the BB has to be chosen 
such that it makes 1 -|- z as small as indicated by (9.5) 
or (9.7), but at the same time the hump is not too wide 
(to avoid large perturbations of the GMB) or too high 
(so that many humps can be in the field of view of the 
observer, to account for the ubiquity of the GRBs). The 
connection between the shape of the hump and the other 
quantities will become clear in Secs. XH D and XV. 

The 5-parameter family of profiles shown in Fig. 1 
emerged in numerical experiments. (For remarks on how 
this shape was arrived at see Appendix G.) Each profile 
consists of two curved arcs and of a straight line segment 
joining them. The upper-left arc (shown in thicker line) 
is a segment of the 4-th degree curve 


r* it - fsf - Apf 


(9.9) 


where tef is given by (8.3). The lower-right arc (also 
shown in thicker line) is a segment of the ellipse 


{r-Bi-Aif {t-tBi-Aof 




An" 


= 1 . 


(9.10) 


The straight line segment passes through the point 
(r,<) = {BiABf + Aq) where the full curves would be 
tangent to each other; the dotted arcs show the parts of 
the curves that are not included in the hump profile. 

The profile in Fig. 1 has five parameters: Aq, Ai, Bq, 
Bi , and xq which determines the slope of the straight seg¬ 
ment. The other quantities are determined by these five. 
Figure 1 is not drawn to scale with respect to the values 





used in actual numerical calculations. In particular, xq 
and Ai are greatly exaggerated."^ 



r -^ 


FIG. 1: Parameters of the bang-time profiles (the profile Is 
drawn not to scale for better readability). See (9.12) for the 
actual values in Model 1 and (9.19) for the values in Model 2. 

This exemplary class of profiles is meant to provide a 
proof of existence of the blueshifting mechanism. Other 
profiles may emerge when modelling the actually ob¬ 
served GRBs; see one example in Subsec. IX C. 

The function B(r) obeys 

2E/r^ ='' -k = 0.4 (9.11) 

and is the same as in the background model, see Sec. 
VIII. A general E would have the form 2E{r) = r‘^[—k + 
J^(r)], where J^(0) = 0, but otherwise is arbitrary [7]. 

One more adjustable parameter defines the path of the 
ray. It is Ate ~ the time-difference between tsiO) and the 
instant when the ray intersects the axis r = 0. The Ate 
determines the coordinate r of the observer who receives 
the ray at t = 0; this way of specifying the initial condi¬ 
tions led to a greater numerical precision. 

A. Model 1 

Three models were constructed. Model 1 is adapted 
to the GRBs of the highest frequency, and satisfies (9.7) 
with some excess. The values of its parameters are 


/Aq\ 


/ 0.00045 NTU \ 

Bq 


0.001 NTU 

Ai 

= 

1 X 10"® 

Bi 


0.03 

\xq J 


\2xl0"i® / 


(the last three parameters are dimensionless). The values 
of Aq and Bq imply the time difference between the max- 


^ The straight segment of the BB is invisible in most of the next 
figures in consequence of xq being extremely small. It was intro¬ 
duced to keep dtg /dr finite between the curved arcs. 


imum of is and its flat part 0.00145 NTU = 1.421 x 10® 
years (~ 0.0103T, with T given by (7.10)). 

The value of Ate in Model 1 that ensured compatibility 
with (9.7) is 

Ate = 1.9949248 x 10"® NTU « 1.955 x 10® y. (9.13) 
This ray reaches the present time t = 0 at 

Hof 

r = Tobsi = 0.76218478306089776. (9.14) 

With (9.12) and (9.13), the blueshift is 

1 + zminb ~ 2.5191965 X 10"® (9.15) 

between the LSH and now. This is slightly better than 
required by (9.7). In this way, the upper limit of GRB 

frequencies in ® Property (1) from Sec. II is accounted 
for. Values oi 1 + z obeying (9.7) could possibly be ob¬ 
tained with an even smaller height of the BB hump (i.e. 
with smaller Aq -\- Bq). However, no precise optimiza¬ 
tion was attempted: this model was meant to be only 
a proof of existence of the blueshifting mechanism ac¬ 
counting for the highest GRB frequencies. More precise 
modelling was done for Model 3, see Sec. IX G. 

The value (9.15) was calculated in two steps, using 
(5.2). First, 1 + z between the LSH and the point Pa of 
coordinates {r,t) = (0,tB(0) + Ate) was calculated; it is 

1 -b = 9.1144891634 x 10"^°. (9.16) 

Then, the redshift between Pa and the observer sitting 
at t = 0 was calculated with the result 

1 -b zoai = 27.63947. (9.17) 

The number in (9.15) is the product of the two. 

The main panel of Fig. 2 shows the radial cross-section 
through the hump in the BB in Model 1, and the ray 
defined by (9.13). The inset shows a closeup view of the 
neighbourhood of tB(0) (top of the hump). The LSH and 
the BB lie so close to each other that they would coincide 
in all the figures, so the LSH is not marked. 

The line that looks vertical in the main panel is part 
of the profiles of the ERH and of the BB that coincide 
only spuriously, because of the scale of the figure. It is so 
close to vertical in consequence of the very small values 
of xq and Ai. Another consequence of the small value of 
Ai is that the lower (ellipse) arc from Fig. 1 is invisible 
here. The dot near the horizontal part of the BB shows 
the point where the ray hits the LSH, see eq. (D2) in 
Appendix D for the exact value of t there (and compare 
(D3) to see how close that point is to the BB). The two 
(seemingly coincident) horizontal strokes above the dot 
show the positions of the ends of the straight segment of 
the BB. The ERH is spherically symmetric around the 
center, so its profile is mirror-symmetric with respect to 
the r = 0 line, but the left part of the profile is omitted 
in the inset, for clarity of the picture. 
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FIG. 2: Main panel: A radial cross-section through the 
BB hump in Model 1. It is of the type shown in Fig. 1, 
with the parameter values given by (9.12). At this scale, the 
ray defined by (9.13) nearly coincides with the BB. See text 
for more explanation. The inset: A closeup view of the 
neighbourhood of the maximum of the BB. The dotted line 
is r = 0. The difference in t between the ray and the BB 
at r = 0, marked by the double arrow, is the value of the 
parameter Ate given by (9.13). 


B. Model 2 


Model 2 was obtained from Model 1 by combining a few 
actions: making the hump lower (by decreasing Aq and 
Bq), changing the slope of its straight part (by increasing 
or decreasing xq), and replacing (9.9) with an analogous 
curve of still higher degree. Each change in Aq or xq had 
to be accompanied by adjusting Ai (to avoid making xi 
larger than Ai, see Fig. 1) and Ate (to keep the ray in 
the range of minimum 1 + z). 

In Model 2, the upper-left arc from Fig. 1 is replaced 
by an arc of the 6-th degree curve, 


{t - tBf - Apf 

Bi^ Bo® 

and its parameters are 


/Aq\ 


/ 0.000026 NTU \ 

Bq 


0.0001 NTU 

Ai 

= 

1 X 10-1° 

Bi 


0.015 

\Xo J 

V 2 X 10-1® / 


(9.18) 


(9.19) 


being the product of two factors analogous to those in 
(9.16) and (9.17). This time they are 

1 + Zais2 = 1.07858890707746014 X 10“^ (9.22) 
l + Zoa2 = 126.246039921. (9.23) 

The 1 -h Zmaxb Is Smaller than the 1 -I- ^max given by 
(9.5), and thereby Model 2 accounts for the lower limit 

of frequencies in ® Property (1) from Sec. II. This ray 
reaches the present time at the radial coordinate 

ritaf 

r = ro 2 = 0.88705643159726955. (9.24) 

The picture corresponding to Fig. 2 is qualitatively sim¬ 
ilar here; it is shown in Fig. 3 to visualise the changes 
caused by replacing the curve (9.9) with the curve (9.18). 



r 


FIG. 3: This figure is the analogue, for Model 2, of the main 
panel in Fig. 2. The meaning of the symbols is the same here. 

Models 1 and 2 are compared in Fig. 4. The BB pro- 
hles consist of the humps from Figs. 2 and 3 surrounding 
the center; further away from the center is constant, 
and so the geometry is Friedmannian. The present time 
is t = 0, and the flat part of ts is at t = tsf given by 
(8.3). The observer lies farther from the hump in Model 
2 than in Model 1. The humps are almost invisible here, 
so they are shown in closeup view in Fig. 5. 

Some details of the numerical calculations for both 
models are described in Appendix D. 


C. Model 3 


Compared to (9.12), Aq is decreased by a factor of ^ 17, 
Bq by the factor 10, Ai by the factor 100 and Bi by the 
factor 2. The height of the hump in the BB is thereby 
decreased from 0.00145 NTU to 0.000126 NTU, i.e. by a 
factor of ~ 11.5. The ray that has the smallest 1 -|- z is 
here defined by 

Ate = 8.83425 x 10“® NTU = 8.657565 x 10® y, (9.20) 
and the blueshift on it is given by 

1 + ^maxb = 1.36167578 x 10"®, (9.21) 


Models 1 and 2 account for the range of frequencies 
of the GRBs, but, as will be seen in Sec. X, do not 
correctly account for their short duration. Model 3 that 
will be presented now is a modification of Model 2 made 
in order to allow us to choose this time interval. 

The modified profile of the BB hump is shown in Fig. 
6 ~ again not to scale. Now the ellipse arc is truncated 
at r = Tge = Bi -I- Xge (ge for gamma-burst end), and 
is replaced by a parabola arc going from r = Cge to r = 
Tge -I- Xag (ag for afterglow) achieved at t = tef given by 
(8.3). In the special case Xag = t/ag = 0 the parabola arc 
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FIG. 4: Comparison of Models 1 and 2. The right-hand half 
of the ERH profile is shown only for Model 1. 

-0.1378 
-0.138 
-0.1382 
-0.1384 
-0.1386 
-0.1388 
-0.139 
-0.1392 
-0.1394 
-0.1396 

FIG. 5: Comparison of the BB profile in Model 1 (upper 
curve, the same as in Fig. 2) with the BB profile in Model 2 
(lower curve, the same as in Fig. 3). 

is absent, and the ellipse goes over (non-differentiably) 
into the straight line t = The parameters i/ge, Xag 
and 2/ag are adjustable, Xge is determined by yge- The 
net effect compared to Fig. I is that the whole hump is 
moved down with respect to the flat part of the BB, by 

toff Ao - yge - 2/ag, (9.25) 



FIG. 6: The bang-time profile of Model 3 and its parameters. 
See text for explanations. 


and the parabola arc interpolates between the ellipse and 
the straight line. The value of Xge will be given in Sec. 
X, while Xag = j/ag = 0 will be assumed in all models. 


X. ACCOUNTING FOR THE LIMITED 
LIFETIMES OF THE GRBS (PROPERTY (2) IN 
SEC. H) 

A. Adapting the model to the observed durations 
of GRBs 

We recall from Sec. II that the observed GRBs typi¬ 
cally last from less than a second to a few minutes [6]. 
Let us take 10 minutes as the reference value. 

If two rays in Model 2 are emitted at r = 0 close to the 
BB, the second one later than the first by 10“^® NTU 
(si 5.15 minutes), then observer O 2 sitting at the ro 2 
given by (9.24) will receive the second ray 1565.56 x 10“^® 
NTU (~ 5.6 days) later than the first. Consequently, 
to get a time-interval of 10 minutes = 1.9414 x 10“^® 
NTU between the received rays at the observer, the time- 
difference between their emission events at the center 
would have to be ^ 1.24 x 10“^® NTU. But rays sep¬ 
arated by such a short interval of t at r = 0, when sent 
back in time, become indistinguishable for a Fortran90 
program at double precision: they (numerically) reach 
the LSH at the same instant with the same 1 -I- 2 :. So, a 
much higher numerical accuracy is needed to trace such 
tiny intervals backward in time. Lacking this, we will 
only show that the duration of a GRB in Model 3 can be 
controlled via the parameter j/ge. 

Numerical experiments with Model 2 showed that 
the program begins to see the difference between past- 
directed rays all the way down to the BB when their 
initial points at r = 0 are separated by 

AT Ri 1.24 X 10"^® NTU 4.4 days). (10.1) 

The earlier of the two rays was chosen such that the 
observer O 2 sees it as (approximately) the first in the 
gamma-ray frequency range. It crossed the LSH at 

<gb2 = -0.13944465218439359 NTU. (10.2) 

The ray that crossed the center later by AT, hit the LSH 
at 

= -0.13944465708554910 NTU < tgb2. (10.3) 

(The ray that crosses the LSH earlier reaches the sym¬ 
metry axis and the observer later, see the figures.) 

For the first set of numerical experiments with Model 
3, called Case I below, we assumed x^g = y^g = 0 and 
took the difference tge 2 — ^Bf as the toff of (9.25). Thus 

Vsesi = ^0 + ^Bf - tge2 = 1.511019508261 X 10“® NTU. 

(10.4) 
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For the ray that defines the beginning of the GR flash 
in Case I, we took the one with 

AG = Qi ='" 8.832869 x 10"® NTU (10.5) 

(see Fig. 2 for the definition of Ate)- With j/ge = j/gesi 
given by (10.4), it reached the LSH at 

tLSH 3 i = -0.13945554226395207 NTU, (10.6) 

and the present time t = 0 at 

r = ro 3 = 0.89000433423592207. (10.7) 

Observer O 3 residing at r = ro 3 will be the new refer¬ 
ence observer. The redshift on the ray defined by (10.5) 


between the LSH and the O 3 world line is 

l + zib = 1.67354 X 10“®, (10.8) 

and satisfies (9.5) with a little excess. 

For the later rays, AG was increased in quanta of 

qt = 10-13 NTU. (10.9) 

The two rays that had AG = and Ate = Qi + 

2qt were still in the gamma-ray frequency range when 
observed, but the one with Ate = qi+ iqt was already in 
the X-ray frequency range, see Table I. 


TABLE I: Limits on duration of GRBs in Model 3 with different values of i/ge and Ate 


Case 

Value of i/ge 

Value of Ate 

Value of T 

1 - 1 - 2 at O 3 

I 

Vgel 

qi -I- 2qt 
qi -I- 3(jt 

9.80519016282838448 x 10"® « 960.90 y 
9.94424450172667973 x 10“® « 974.54 y 

1.67815 X 10“® 
5.05133912545 x lO"® 

II 

Vgell 

qi + 2qt 

qi + 

9.86541227093332829 x 10"® « 966.81 y 
1.00044488426559296 x 10“® « 980.44 y 

1.676761 X 10"® 
5.051339211 X lO”® 

III 

Vgelll 

qi + 2qt 

qi + 

9.92591136595910132 x 10“® « 972.74 y 
1.00649385349405952 x 10“® « 986.36 y 

1.67476 X 10“® 
5.051339288 x lO"® 


In units of NTU, in Case II the value of is 

2/ge3ii = 2/ge3i + 2^* = 1.511019528261 x lO"®, (10.10) 
and in Case III it is 

2/ge3iii = 2/ge3i + 4<7t = 1.5II0I954826I x lO"®. (10.11) 

The values of these and other parameters for all cases 
are given in Table 1. Assuming that all rays in the table 
were emitted in the visible range, in each case the ray 
with 1 -|- 2 ~ 10 “® is observed in the gamma-ray range, 
while the one with 1 -|- z ~ IQ-^ is observed in the X-ray 
range [25], i.e. after the end of the CRB. The parameter 
T is the t-coordinate of the instant when O 3 sees the 
ray. Values of T are given both in NTU (the first one) 
and in y. The first line of 'T in each case provides the 
lower limit, the second one provides the upper limit on 
the duration of the CRB. By this criterion, Table I shows 
that the duration of the CRB in Model 3 increases when 
?/ge increases, so it is controlled by j/ge. (In Model 2, 
which is the limit j/ge = Aq of Model 3, the gamma-ray 
“burst” would last more than 10 ^ years.) 

The values in Table I are only rough estimates of the 
CRB duration because they assume that in all three cases 


the CRB begins at t = 0, while in fact the observed 
frequency might have entered the gamma range earlier. 
Without a precise identification of the onset of a CRB in 
a model, no better estimate is possible. 

The estimates given here ignore one more problem that 
is important for detecting the GRBs, namely their inten¬ 
sity. For a CRB to generate a signal in a detector, it must 
not only be in the gamma range, but also must be suffi¬ 
ciently strong. The problem of calculating the intensity 
of a CRB is left for future consideration. 


B. The time-profile of observed frequencies in case 
I of Model 3 

We will now follow the sequence of radial rays received 
by observer O 3 in case I of Model 3, from the (approxi¬ 
mately estimated) moment when the BB hump first ap¬ 
pears in her field of view up to the moment (also approxi¬ 
mate) when it disappears. The parameters of several rays 
emitted at characteristic points of the LSH are given in 
Table H. The values of t are given in NTU unless speci¬ 
fied otherwise. 
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Selected rays are shown in Fig. 7 along their whole 
paths, and in Fig. 8 in the neigbourhood of the BB. 
Note how the rays abruptly change their slopes when 
crossing the steep wall of the ERH but remain smooth 
on intersection with the other branch of the ERH. This 
is because the near-vertical part of the BB hump creates 


a quick rise in the ERH from t = 0 to a large value above 
the upper margin of the figure. This is not a general 
property of all ERH, but only of this particular model. 

The spurious intersections of different rays in Fig. 8 
are caused by a limited resolution of the figure; in reality 
the rays do not intersect. 


TABLE II: Parameters of the rays shown in Figs. 7 and 8 


Ray label 

(t, r) at the LSH 

t at O 3 

1 -|- 2 : at O 3 

Rl 

-0.13944546999999999 

0.0150000000139671984 

-0.02621663 NTU 
ss 2.569 X 10^ y ago 

445.36 

R2 

-0.13941000000000001 

0.0149998605746265272 

-0.01342874 NTU 
ss 1.316 X 10® y ago 

130.439 

R3 

-0.13935 

0.0142870667194671501 

-3.4317801 X 10”® NTU 
« 3.363 X 10® y ago 

75.9982 

R4 

-0.13933632431674814 

0.0075 

-1.054319 X 10"® NTU 
« 1.033 X 10® y ago 

92.142 

R5 

-0.13933556296232263 

0.0 

-5.50955769 x 10"^ NTU 
« 5.399 X lO’’ y ago 

143.9397 

R6 

-0.13933607609165491 

0.00700320078026953671 

-6.484649 x 10"® NTU 
« 6.3549559 x 10® y ago 

203.5369 

R7 

-0.13937350162416995 

0.0147673852333973319 

-6.113304 x 10"® NTU 
« 5.991 X 10® y ago 

1.6002 

R8 

-0.13939833897747514 

0.0149860617889288825 

« 2 X 10"® NTU (*) 

« 1960 y ago 

0.097862 

R9 

-0.13945554226395207 

0.0150000000186084614 

0 (now) 

1.67354 X 10"® 


(*) - see remarks in text. 


Ray RI is emitted at the LSH at nearly the same r that 
is marked with a dot in Fig. 3. This is one of the earliest 
signals from the neighbourhood of the hump that crosses 
the world line of observer O 3 . At a short initial segment 
of its path this ray acquires some blueshift. However, 
on exit from the ERH the accumulated blueshift is not 
sufficiently strong to survive and is completely obliter¬ 
ated by later-acquired redshifts. In the end, as seen in 
Table H, ray Rl reaches observer O 3 with a redshift suf¬ 
ficient to shift the lowest-frequency visible light to the 
high-frequency segment of microwaves, and the highest- 
frequency visible light to infrared. 

The 1 -I- z at O 3 decreases from RI to R3. This is 
because, when the emission point of the ray is moved up 
the hump, the ray travels a longer distance inside the 
ERH, so it acquires more blueshift. Moreover, the travel 
time of the ray outside the ERH becomes thereby shorter, 
so the accumulated redshift becomes smaller. 


On rays received after R3, the observer sees increasing 
1 + z until long past ray R5, emitted at the center of 
the hump. This happens so because somewhere between 
R3 and R4 the rays begin to stay within the ERH for 
shorter and shorter times. At r = 0, where R5 has its 
origin, the ERH is tangent to the BB. Consequently, the 
initial point of R5 is outside the ERH and the ray builds 
up positive redshift from the very start.® On rays later 
than R5, I -f z at O 3 still increases for a while. This is 
because with increasing difference in r between the light 
source and the observer, the redshift acquired outside the 
ERH outbalances the blueshift acquired inside the ERH. 


The redshift decreases locally on all rays when they pass through 
the wall of ERH, seen in Fig. 4. But this has little effect on the 
final value of 1 -\- z. 
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FIG. 7: Selected rays from Table II. The rays between R3 
and R9 are not shown for lack of space. The arrow points 
to the center of the BB hump. The upper horizontal straight 
line is the present time. See text for more details. 



- r r -^ 

FIG. 8: Closeup view of the neighbourhood of the BB hump 
in Fig. 7. Initial points of the rays listed in Table II are 
marked with dots. The vertical straight line is the axis of 
symmetry. For rays between R4 and R9 only their initial 
points are marked, for lack of space. 


Somewhere between R 6 and R7 1 + z at O 3 achieves 
a maximum and then begins to go down. Between these 
two rays there is a discontinuity that makes it extremely 
difficult (perhaps impossible) to hit, with past-directed 
rays originating at r = 0 , the middle of the arc between 
the points marked R 6 and R 8 . The rays hit the LSH ei¬ 
ther above that middle point, with 1 -|- z ~ 1.5 at O 3 , or 
below it, with 1 -|-z « 0.02 at O 3 . This is where observer 
O 3 begins to see blueshifts. If this is not a numerical ef¬ 
fect, but a reflection of a real discontinuity in the model, 
then blueshifts appear abruptly, and are initially moder¬ 
ate: 1 + z = 0.02 takes the visible range into the UV. 

A similar instability exists between rays R 8 and R9: 
the change in 1 -|- z measured by O 3 from the order of 
10“^ to 10“® occurs abruptly. If both rays are emitted 
as visible light, then ray R 8 is observed still in the UV 
range, wile R9 is already in the gamma range. Ray R9 
is received by observer O 3 at t = 0, i.e., now. 

Again, if this is a reflection of a real discontinuity in the 
model, then the gamma-ray flash would appear suddenly. 


without a continuous transition from the UV range. This 
seems to agree with the observations of the real GRBs [ 6 ]. 

The asterisk (*) in Table II, in row R 8 and column “t 
at O 3 ”, indicates that the t at O 3 for Ray 8 had to be 
hand-corrected. The time difference between rays R 8 and 
R9 at O 3 is so small that numerical roundoff errors inter¬ 
fere with it. The numerically integrated ray R 8 overshot 
the ro 3 given by (10.7), and the final t reported by the 
program was greater than zero, which is impossible (if 
this were true, then R 8 and R9 would have to intersect 
somewhere). So, the t corresponding to r = ro 3 on ray 
R 8 had to be estimated by comparing the differences in 
t between R 8 and R9 at other values of r. Near O 3 they 
were ~ 2 x 10"® NTU. 

For rays received by O 3 later than R9, the familiar in¬ 
stability shows up once more:® the transition from the 
gamma-ray range to the X-ray range occurs abruptly 
when At is increased by 10“^® NTU; see Sec. XA. This 
is the transition from the GRB to the afterglow. 

In addition to the abovementioned discontinuities that 
are justified by the properties of the BB profile, there 
is one more instability that seems to be of purely nu¬ 
merical character, and which this author was not able to 
explain or remove. Namely, when the rays from Fig. 8 
are retraced back in time from the observer, they coincide 
satisfactorily with the original rays only up to the ERH. 
A numerical instability at the very high and steep wall of 
ERH causes that the change of position of the end point 
(in the past) of the ray is a discontinuous function of the 
initial position. The (so far irremovable) discrepancies 
between the initial point of the future-directed ray and 
the end point of the past-directed ray were of the order 
of 10"® NTU at r = 0. 


XI. ACCOUNTING FOR THE AFTERGLOWS 
(PROPERTY (3) IN SEC. II) 

For this section, we use Case I of Model 3. 

In this model, the afterglow appears necessarily. The 
first ray in the afterglow is the one listed in the second line 
of Table I. (Figure 9 does not show it; it would mostly 
coincide with the BB at this scale.) If it was emitted in 
the visible range, then it is observed in the X-ray range. 

As time goes by, the observed frequency goes down. 
For example, the earlier of the two rays shown in Fig. 9 
reaches observer O 3 with z = —0.0008107, thus its ob¬ 
served frequency is nearly the same as the emitted one. 
The later ray reaches O 3 with z = 0.598, which shifts 
nearly the whole visible range into the infrared. The 
points marked “1” and “3” in the figure are at the lo¬ 
cations where the earlier ray hits the LSH and the BB, 


® The discontinuities in 1 -|- z might be caused by the non¬ 
differentiability of the BB profile at r = B\ — xg, r = Bi + x\ 
and r = Bi + Xge', see Figs. 1 and 6. 
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respectively, the points “2” and “4” show the correspond¬ 
ing events for the later ray. However, this frequency drift 
occurs at the cosmological time-scale and would not be 
observable with the current technology [26]. So, it cannot 
be responsible for the gradual fading of the afterglow. 



FIG. 9: Two characteristic rays in the afterglow sector of 
Model 3. See explanation in the text. 


The factor that is responsible for the fading of the sig¬ 
nal is the intensity of the received radiation. When it 
drops below the sensitivity of the detector, the afterglow 
is, in the technical sense, ended. The question of calcu¬ 
lating the intensities of the afterglows in this model is 
left for a separate investigation. 

Model 3 does not account for the duration of an after¬ 
glow quantitatively, and does not contain any parameter 
that could control this period. As can be seen from Table 
I, observer O 3 in Case I would begin to see the afterglow 
at t ^ 974.54 y after now. The later ray in Fig. 9 is still 
well within the afterglow, and reaches O 3 at 3.5714 x 10“® 
NTU ~ 3.5 X 10® y after now. So, assuming that the in¬ 
tensity of the ray would still be sufficient for detection, 
this model-afterglow would remain visible for observer 
O 3 for nearly 350 000 years, while the longest-lasting af¬ 
terglows observed in reality are visible only for several 
hundred days [27]. Thus ® Property (3) from Sec. II is 
accounted for only qualitatively. 

The blueshifts in the afterglow of Model 3 arise when 
the rays pass through the ERH wall seen in Fig. 9. Con¬ 
sequently, in order to reduce the duration of the afterglow 
one should force the ERH to stay nearer to the BB hump, 
so that the observed rays begin to bypass the ERH as 
early as possible. This should be possible with BB pro¬ 
files having more parameters, but is not possible in Model 
3. Other BB profiles would also be needed to account for 
lower-frequency afterglows. With those profiles, the path 
of the ray inside the ERH should be suitably short so that 
the ray would acquire less blueshift. 


XII. NONRADIAL RAYS 

A. An exemplary nonradial ray hitting the BB 

As an illustration to the remark under (4.8), Fig. 10 
shows the behaviour of a nonradial ray that hits the BB 
where dt^/dr ^ 0. This is a projection of the ray and 
of the edge of the BB hump from Fig. 8 on a surface 
of constant t along the flow lines of matter (in comoving 
coordinates the image is the same at every t). For brevity, 
from now on we will use just the word “projection” to 
denote this kind of image. The coordinates in the figure 
are related to (r, (p) by 

X = rcosp, y = rsimp. (12.1) 

The outer edge of the hump is marked with the big cir¬ 
cle, its center with the cross. The ray is emitted to¬ 
ward the past at a point slightly off-center (r = 0 . 0001 ), 
in a direction tangent to r = constant, at the same 
t = —0.13933160382638388 NTU as ray R9. It meets 
the BB at the steep slope of the hump. As predicted, it 
becomes nearly tangent to an r = constant surface before 
hitting the BB. The dot ^ 1 marks the point where the 
ray crosses the LSH (with 1 + z = 0.71003236199623432 
relative to the initial point at r = 0), the dot # 2 marks 
the point where it crosses the BB. Somewhere between 
these two points the calculation becomes unreliable be¬ 
cause of numerical errors (see Appendix E). The inset is 
a closeup view on the initial point of the ray; it shows 
how far off the center the ray begins. 



FIG. 10: A nonradial past-directed ray that aims at the BB 
where dts/dr 7 ^ 0. It becomes tangent to an r = constant 
surface at the BB. See text for more explanation. 
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B. Angular size vs. C 

The constant C that first appeared in (4.7) is a measure 
of non-radialness of a ray. It can be related to the angular 
size of an object grazed by the ray. It follows from (4.13) 
using (5.3) that at the point where dr/dA = fc’’ = 0 
(i.e. where the ray becomes tangent to an r = constant 

sphere), C, R Rq and z zq obey 

i?o = |C|/(I + zo)- (12.2) 

The i?o determines the impact parameter, see Fig. 11. 



FIG. 11: Calculating the angular size of an object. See text 
for explanation. 

Figure 11 shows the projection of a ray, similar to Fig. 
10. O is the position of the observer, D is the center of 
the BB hump, the dotted circle is the edge of the hump 
and the smaller solid circle represents the object whose 
angular size we wish to calculate. The point P is the 
event where (12.2) holds, P’ is the intersection of the solid 
circle with the circle that has the center at O and radius 
OD, and B is where the (past-directed) radial ray OD 
would enter the L-T region. This is an illustration only; 
inside the L-T region the projection of a real nonradial 
ray would not be straight - see next pictures. 

To calculate the angular radius a of the solid circle 
in Fig. 11 as seen by the observer at O, two approxi¬ 
mations must be made. In a flat space, a would be the 
length of the arc DP’ (call it Ia) divided by the length of 
the straight segment OD (call it ij^). In the cosmological 
context, for in we can take the angular diameter distance 
(ADD). In the L-T model (3.1), for an observer at the 
center of symmetry, ADD is the value of R taken at the 
location (t,r) of the observed object. In the Friedmann 
limit (7.1), R becomes rS{t)^ and every point is a center. 
The difficulty in the present case is that the spacetime is 
not spherically symmetric around O, and there is no op¬ 
erational definition of the angular diameter distance in a 
general cosmological model. So, the two approximations 
to make are these: 

1. We take Rq for Ia- In a flat space, the difference 
between Rq and Ia would be £/j(a — sina), so the relative 
error at oq = 1 degree = 0.01745 rad (the actual angular 
radius of the BB hump for O, see below) would be I — 
(sinao)/ao ~ 5 x 10“^. 

2. To calculate in we assume that the ray OD prop¬ 
agates in the Friedmann background all the way from 


r = ro 3 to r = 0. The error induced by this assumption 
can be estimated in two ways: 

(a) As measured by the values of the r-coordinate the 
segment BD is ~ 0.017 of OB. 

(b) The center r = 0 is reached by a ray in Model 3 at 

= -0.13932806655865945 NTU, (12.3) 

and in the Friedmann model at 

tp = -0.13934135010087942 NTU, (12.4) 

so the difference is 

tM3 - tp = 0.00001328354221997 NTU, (12.5) 

which is ~ 0.0095% of |tM 3 |- 

Under assumptions 1 and 2, the ADD from O to D is 

in = 3.45912604264922777 x 10“^ NLU; (12.6) 

see Appendix F for derivation. The angular radii in col¬ 
umn 3 of Table III were calculated as Ro/in, where Rq 
was a by-product in computing the path of the ray. 

C. Nonradial rays received by observer 03 in 
Model 3 

Now we will display several nonradial rays received by 
observer 03 at the present time from different directions. 
We will follow them back in time from the initial point 
at 03. Table III lists their parameters. The angular 
radius is calculated as explained above. In Figs. 12 
to 15 the horizontal coordinate a; is r cos ip and the la¬ 
bels of the rays are the same as in Table III. Note that 
the redshifts at the LSH increase when the impact pa¬ 
rameter increases. The radial ray would have the strong 
blueshift given by (10.8); the changeover to redshift oc¬ 
curs abruptly as soon as the ray ceases to be radial. Ray 
6 has larger redshift than the background (8.5); this is 
probably a peculiarity of Model 3 caused by the steep 
rise of the BB hump profile. 

Figure 12 shows the projection of the rays from Table 
III in a vicinity of their end points. Fig. 13 is a closeup 
view on the L-T region over the BB hump. The labels of 
the rays are printed at their end points; at the left margin 
of both figures rays 1-8 are ordered from bottom to top. 
The dotted circle is the edge of the BB hump, the cross 
marks the center r = 0. The large dots in Fig. 12 mark 
the points where the rays intersect the LSH. The end 
points of the rays are where the numerical calculation 
determined that the ray crossed the BB. 

The rays abruptly change their direction every time 
they intersect the steep wall of the ERH (see Figs. 2 
and 3). The change is sharper when the ray is closer to 
the BB; this happens on the final exit from the ERH. 
The gentle deflections take place when the rays travel 
between the two branches of the ERH. (The ERH cannot 
be shown in these figures: in such a projection the steep 
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wall would coincide with the hump circle, and the other 
branch projects onto the whole disk inside the circle.) 

The angle of deflection of a ray depends on the interval 
of t that the ray spends between the two branches of the 
ERH. Ray 1 meets the ERH nearly head-on and does 
not strongly change direction on first encounter (after 
the first encounter with the ERH ray 1 stays between 
the two branches only briefly, similarly to ray R9 in Fig. 
8). However, after the second encounter, it stays between 
the two branches for a longer interval of r (cp. Fig. 3). 
Then it is close to the BB, where R is small, and is forced 
to bend around in agreement with (4.8). 



FIG. 12: View on the rays listed in Table III as they fly over 
the BB hump. See explanation in the text. 


The other rays meet the ERH at smaller angles than 
ray 1, so they stay between the ERH branches for longer 
times. For rays 2 and 3 this gives the effect that they are 
deflected much stronger than ray 1. For ray 4 and the 
other ones, a different effect prevails: they stay farther 
from the BB, so approach the BB later and stay within 
its influence for a shorter time, therefore the deflection 
angle decreases again. Rays 7 and 8, which do not enter 
the L-T region, do not feel its influence and propagate 
just as in a pure Friedmann spacetime. 

Figure 14 shows the projections of the rays from Fig. 
12 on the plane y = 0. The intersections of the BB and of 
the ERH with this plane are also shown (the right part of 
the ERH profile is suppressed to avoid clogging the im¬ 
age). In this projection, rays 1 to 4 very nearly coincide 
from the left margin of the figure up to the neighbour¬ 
hood of the right shoulder of the BB hump, and then 
they split in consequence of the different deflection an¬ 
gles. Rays 7 and 8 nearly coincide all the way. Ray 6 is 
clearly visible only between the “knot” where the projec¬ 
tions of all the rays intersect and the right wall of the BB 



FIG. 13: A closeup view of the region around r = 0 in Fig. 
12 . 


hump, elsewhere it nearly coincides with 7 and 8. Ray 5 
is visible on both sides of the knot, and to the right of 
the BB hump, where it stays close to 7 and 8. 



X -s- 


FIG. 14: Projections of the rays from Fig. 12 on the radial 
plane y = 0. See explanation in the text. 


D. Relation between the GRBs and the CMB 

Redshifts at the LSH listed in Table HI must be com¬ 
pared with the redshift in the background model. The 
result (8.5) was obtained by substituting numerical val¬ 
ues of model parameters into exact formulae. Numerical 
integration of the null geodesic equation gives, for ray 7, 
a result differing by ~ 0.07%: 

I -p = 951.91469714961829 (12.7) 
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TABLE III: Parameters of nonradial rays 


Ray 

C 

ang. radius (°) 

r at LSH 

1 + z a.t LSH 

r at BB 

1 

0.0001 

0.0126947 

5.06068826958865417 x 10“^ 

287.16608259998554 

7.95216518223493263 x 10"^ 

2 

0.0005 

0.064259 

4.85198622293232865 x 10~^ 

334.76688994090046 

7.93592114807387117 x 10"^ 

3 

0.001 

0.1337447 

4.52659901407877416 x 10“^ 

460.87136026281968 

7.45576021007735151 x 10"^ 

4 

0.0016 

0.226233 

4.79104086556434161 x 10“^ 

703.10147012965876 

7.66030728333577382 x 10"^ 

5 

0.0032 

0.4743 

5.45954223781688272 x 10“^ 

945.43487410592388 

8.45197957881479167 x 10"^ 

6 

0.0045 

0.6726575 

5.76620501069033120 x 10~^ 

970.10190933005345 

8.60372698978428968 x 10"^ 

7 

0.007 

1.00097 

5.95657733949151447 x 10“^ 

951.91469714961829 

8.96185695386919057 x 10"= 

8 

0.008 

1.1452 

6.02258192362461128 x 10“^ 

951.91132098857997 

8.86836611484453224 x 10"^ 


(see Table III for ray 8 ). Improving the consistency be¬ 
tween (12.7) and (8.5) would require a higher numerical 
accuracy, and the integration time would become pro¬ 
hibitively long. The value (12.7) was found with the 
step in the affine parameter being AA = 10“® up to 
X = Xb '= — 0.03, and 10“^^ for x > Xb- 

As Figs. 13 and 14 show, ray 7 nearly grazes the BB 
hump. Column 3 in Table III then implies that the area 
of the sky with redshifts different from the background 
would hll a cone of angular radius ^1° for observer O 3 . 
This is twice the resolution of current GRB detectors. 
Thus, rays 1 to 5 would be hidden within this circle. 
Checking for their presence would require aiming a de¬ 
tector at this area while the GRB is still on. 

With the redshifts given in column 5 of Table III, all 
those rays, if emitted in the visible range, would be seen 
by the observer in the microwave range [25]. It must be 
recalled here that in Model 3 (and also in Models 1 and 
2 ) blueshifts visible to the observer are generated only 
on radial rays, so the gamma-ray signal, including the 
afterglow, would have a strictly point source. 

The only reason to worry are rays between 5 and 7, 
on which the redshift is different from that of the CMB, 
and which would be visible around the GRB signal in 
the present version of Model 3. This shows that Model 3 
would need a modification of its parameters (which is nec¬ 
essary also for other reasons) to hide the lower-redshift 
signal within the unresolved patch of gamma radiation. 
With that improvement done, the model will not pre¬ 
dict perturbations of the observed CMB larger than the 
GRBs actually cause. If more-exact measurements in 
the future detect some variability within the presently 
unresolved gamma-ray dots, then it will be possible to 
compare it with the model prediction and improve the 
model accordingly (or discard it). 


The sources of the GRBs are seen as fuzzy circles of about 1 
degree in diameter, private communication from Linda Sparke. 


XIII. DEALING WITH THE COLLIMATION OF 
THE GRBS (PROPERTY (4) IN SEC. H) 


The blueshifts that account for the GRB frequency 
range occur on radial rays, and these are emitted isotrop¬ 
ically by the L-T region. Thus, there is no real collima- 
tion of the GRBs in our models. However, an observer 
may have an illusion that they are collimated. Namely, 
as seen from Table III, rays reaching the observer at an¬ 
gle /3 away from radial are seen with large redshifts down 
to less than /3 = 0.0127°. So, the GRB source appears to 
the observer as nearly point-like. 

To account for anisotropy and possibly for the colli- 
mation, one would have to use an anisotropic model for 
the BB hump.^^ A good candidate is the quasi-spherical 
Szekeres (QSS) metric [28, 29], [7]. It contains the L- 
T model as a subcase, but in general has no symmetry. 
Nevertheless, it can be matched to any Friedmann back¬ 
ground just as well. 

Blueshifts visible to the observer are generated in an 
L-T region only on radial rays, but there is no obvious 
definition of a radial direction in a QSS spacetime. How¬ 
ever, this spacetime contains a flow line of matter that is 
a natural generalisation of the center of symmetry exist¬ 
ing in the L-T models; it was called origin in Ref. [30]. 
So, light rays passing through the origin in a QSS space- 
time should share some properties with the radial rays of 
the L-T models. Gonsequently, it has to be verified what 
happens with blueshifts when a QSS model is employed 
(a problem that deserves an investigation independently 
of the GRB context), and then the question of collima- 
tion of the GRBs can be reconsidered. 


But the collimation of the real GRBs is a hypothesis, not an 
observationally verified fact [6]. 
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XIV. ACCOUNTING FOR THE LARGE show details of the main panel in two regions where the 

DISTANCES TO THE GRB SOURCES curves form particularly complicated knots. 

(PROPERTY (5) IN SEC. H) 


The distances to the GRB sources are inferred by mea¬ 
suring the redshifts in the afterglows [31] and assum¬ 
ing that the redshift-distance relation that holds in the 
Robertson - Walker models applies to them. However, 
blueshifting renders the redshift-distance relations multi¬ 
valued (see below and also Ref. [20]), so when blueshifts 
are present along the ray, redshift fails to be an indicator 
of distance. 

In the models hitherto presented, the sources of the 
GRBs lie close to the Big Bang, so they are ^ 13 x 10® 
years to the past from the present observers - still farther 
than the redshift measurements imply, and in this way 
® Property (5) is accounted for. 

Figure 15 shows the relation between the redshift 2 ; 
along the rays from Table III measured by observer O 3 
(who is placed at a; = —roa given by (10.7)) and the coor¬ 
dinate X = r cos (p of the light source. The different values 
of z at which the curves have their end points result from 
different levels of numerical “approximation to infinity” 
- on every nonradial ray the redshift between the BB and 
any observer must be infinite. The background redshift 
at the LSH given by (12.7) is marked with the horizontal 
line, the large dots mark the value of z between the LSH 
and O 3 on each ray. The redshift profiles on rays 6 , 7 
and 8 coincide at the scale of the figure, except for the 
end points. Profile 6 ends at (r, z) ~ (0.086,122 303.73), 
profiles 7 and 8 end at (r, z) « (0.0896,299 116.45) and 
(r, z) ~ (0.0887,166 586.96), respectively. The graph be¬ 
gins at X = —0.5 because at a; < —0.5 all the curves very 
nearly coincide and behave in the standard way: z just 
increases along each ray. 



FIG. 15: The relation between the redshift along the rays 
from Table III registered by observer O 3 and the coordinate 
X = r cos ifi of the light source. See explanation in the text. 


The main panel in Fig. 16 shows a closeup view on 
the neighbourhoood of (x, z) = (0,100) in Fig. 15, where 
the redshift behaves in complicated and untypical ways. 
The dotted vertical lines mark the values of x at which 
the outer ERH intersects the plane y = 0- The two insets 




FIG. 16: Gloseup views on the region around (x, z) = (0,100) 
in Fig. 15. The redshift is seen to be non-monotonic along 
the rays inside the ERH. See explanation in the text. 


In order to understand these graphs, one has to read 
them in the correct way. The observer is at x = —ro 3 
given by (10.7), beyond the left margin. One has to fol¬ 
low the curves beginning at the position of the observer 
in the direction of increasing affine parameter. As long 
as the rays are away from the BB hump (which has its 
edges at X ~ ±0.015), all curves very nearly coincide 
and z increases along them. When a ray goes more than 
halfway around the hump center (cp. Fig. 13), x begins 
to decrease. This is why each of the curves 2, 3 and 4 
turns back at a certain point (and so does curve 1 near 
its final exit from the ERH, but inconspicuously). 

Let us follow ray 2 as an example. Redshift along it 
increases until the ray crosses the outer branch of the 
ERH for the first time. At that crossing, z has a local 
maximum (z ~ 84, see left inset) and begins to decrease. 
It decreases until the ray crosses the inner branch of the 
ERH for the first time (at x ~ —0.014 with z ~ 79, left 
inset). Then z goes up until the ray crosses the inner 
ERH for the second time (at x ~ 0.0118 with z Ri 212, 
right inset; the graph should be viewed in colour and 
magnihed). Erom here, z decreases down to the minimum 
achieved when the ray crosses the outer ERH for the 
second and last time (at x « 0.007 with z « 80, the 
main panel). Beyond that point, z keeps increasing until 
the ray crosses the LSH; the coordinates of that point are 
given in Table HI; see also Eig. 15. 

Two local maxima and two local minima can be seen 
along each of rays 1-5. Ray 6 does not intersect the 
ERH, so has no local extrema, but redshift along it shows 
slight departures from the background profile. Rays 7 
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and 8 that propagate in the Friedmann background do 
not feel the presence of the ERH. 




FIG. 17: Upper panel: The (z,r) relation along the rays 
from Table III. The observer is at r = ro 3 beyond the right 
margin. Note the radial ray labelled “0” - the only one that 
displays blneshift to the observer. The inset shows that 2 : = 
— 1 at the BB. See more explanation in the text. Lower 
panel: A closenp view of the knot in the npper panel. 


Figure 17 shows the {z,r) relation in a neighbourhood 

def 

of the value of r = terh, at which the rays intersect 
the outer branch of the ERH. The observer is at r = roa 
beyond the right margin. This figure has to be read in 
the same way as Fig. 16, i.e., beginning at the observer 
and going in the direction of increasing affine parame¬ 
ter. Outside the ERH, z increases with decreasing r for 
1^03 > r > terh, and increases with increasing r when r 
again becomes greater than terh , going to infinity at the 
intersection of the ray with the BB. The radial ray (la¬ 
belled “0”) is also included; it is the only one on which 
the observer sees z < 0 for a source at the LSH. The 
small inset in the main panel shows z(r) on the radial ray 
near the BB, where z —— 1; the vertical stroke marks 
^ = ?'ERH- The lower panel shows a further-magnified 
closeup on the near vicinity of r = terh- 

These graphs demonstrate that the ERH, which was 
determined using only radial rays (see Sec. VI and Ref. 
[5]), is the locus of extrema of z also on nonradial rays. 
Finally, Fig. 18 shows an example of a deception that 


may befall an observer when she uses a Friedmann model 
to interpret the redshift that was generated in an inho¬ 
mogeneous Universe. This is the z(r) graph along the 
later of the two (radial) rays shown in Fig. 9. The r- 
coordinate of observer O 3 is marked with the right dot 
in the main panel. The r-coordinate of the LSH is marked 
with the left dot. The redshift measured by O 3 for light 
emitted at the LSH is z^ « 0.598. When interpreted 
against the standard ACDM model with the parameters 
given in (7.8), this leads to the conclusion that the ray 
was emitted 5.9 x 10® years ago [32]. The observer would 
not suspect that between the LSH and herself (i.e., closer 
than the dot in the graph) there are sources of light for 
which the redshift would be much larger, up to ^ 230. In 
our model this ray was emitted at the LSH 0.13945 NTU 
= 1.3666 X 10^® years ago. The redshift for a ray emitted 
that early in the ACDM model would be between 25 and 
26 according to the same “cosmology calculator” [32]. 



FIG. 18: An illustration to the deception that may happen 
when redshift is interpreted against a Friedmann background 
model. See text for explanation. Main panel: The full 
redshift profile between the BB (left end of the curve) and 
the observer (right dot). The left dot marks the redshift at 
the LSH. Inset: A magnified view of the neighbourhood of 
the left dot in the main panel. 


XV. ACCOUNTING FOR THE MULTITUDE OF 
THE GRBS 

So far, we discussed models of a single CRB. A model 
that would account for the multitude of observed GRBs 
should be imagined as a Friedmann background contain¬ 
ing many humps like the ones in Figs. 2 or 3 of different 
shapes, spatial extents and heights above the flat part of 
tBir), placed at different comoving positions. 

For estimating how many GRB sources in the sky our 
model could accommodate, let us assume that all sources 
are described by Model 3. As explained in Sec. XH, the 
GRB generator in Model 3 would fill the angular diame¬ 
ter of 2 degrees for observer O 3 . Consider a unit sphere 
and a cone with the vertex at the center of the sphere 
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that has the opening angle 2i9o. This cone subtends the 
solid angle 27r(l — cos-do). The full solid angle is iw. So, 
counting naively, with do = 1° = 0.01745 radians, the 
number of the corresponding solid angles in the full solid 
angle would be 

A/'-y = 2/(1 - cosdo) = 13 131. (15.1) 

However, to get a more realistic estimate we have to con¬ 
sider each circle on the unit sphere being inscribed into a 
square, and then calculate the ratio of the full solid angle 
to the area of that square. On a plane, the surface area 
ratio between a square of edge a and a circle of diameter 
a is 4/7r. So, the number above must be multiplied by 
7r/4, and then the (approximately estimated) number of 
GRBs that could be visible in the sky at the same time 
comes out to be 

= 10 313. (15.2) 

This is probably too few to be realistic. However, with 
the angular size of the GRB generator reduced by the 
factor /, the number A// would be multiplied by /^, and 
Model 3 needs such an improvement anyway. 

XVI. POSSIBLE WAYS OF IMPROVING THE 
MODEL 

The influences on 1 -|- 2 of the various parameters of 
Models 1, 2 and 3 are interconnected. For example, de¬ 
creasing xq alone (see Fig. 1) has the immediate result of 
mcreasing 1 + z (because the past-directed ray then hits 
the BB behind the hump and displays larger z to the ob¬ 
server). A decrease in 1 -I- ^ is achieved when this change 
in xq is accompanied by decreasing Ai (to make the steep 
straight segment longer) and decreasing Ate (to make 
the ray hit the hump where it is steep, and at t being 
as early as possible). But if the straight segment is too 
steep, then the ray and the LSH become nearly parallel 
and the numerical program has difficulty locating their 
intersection. In consequence of such inter-dependences 
between parameters, decreasing the height and width of 
the BB hump while keeping 1 + z sufficiently small is an 
extremely tedious process that requires a great number 
of long-lasting numerical experiments. (This is why the 
background model was left in an imperfect form, in which 
the value of z at the LSH, eq. (8.5), does not agree with 
(7.11). This point needs to be improved, too.) Thus, the 
whole optimization should best be done by a computer 
program. Optimization by hand, applied in this paper, 
already allowed for a radical decrease in the size of the 
BB hump,^^ but it is not known how much smaller the 
hump can still be made. 


The crucial parameters are Aq, Bq and Bi that determine the 
height and width of the hump; Ai is also important, but it is 
already small, see (9.19). Obtained by the hand-optimization 


With a lower BB hump, also the duration of the GRB 
would be reduced, and the hump would be further away 
from the observer, automatically making the angular size 
of the GRB generator smaller. With narrower humps, 
more of them could be fitted into the full solid angle. 

One way of improving the model is to put more param¬ 
eters into the BB profile. However, some of the improve¬ 
ments were achieved by nontrivial means. For example, 
surprisingly good results were achieved by taking higher- 
degree curves instead of ellipses in the profile shown in 
Fig. 1, see Appendix C for more on this. Computerized 
optimization could help also here. 

Higher numerical accuracy will be needed for control¬ 
ling the duration of the GRBs and of the afterglows in the 
model. The discontinuities reported in Sec. X are clearly 
related to insufficient accuracy. However, a higher accu¬ 
racy will result in longer run-times of the programs, and 
this is one more reason why optimization by computer 
would be useful. 

The observed durations of the GRBs are determined 
by the intensities of the gamma-ray flashes. The ways of 
calculating the intensities must still be found. 

It would be desirable to have a hump profile whose 
ERH would not include a wall as high and wide as that 
in Models 1 ~ 3; see Figs. 2 - 4, 8 - 9 and 14. As 
can be seen from (6.2) and (6.3), t on the ERH is large 
where \dtB/dr\ is large. Thus, the ERH wall is narrower 
when the segment of large \dtB/dr\ is shorter, and is 
lower when the value of \dtB/dr\ is smaller. A lower 
and thinner wall would result in reducing the duration of 
the GRB, and could help in controlling the duration of 
the afterglow, both of which are too long-lasting in the 
current models. However, a smaller width of the ERH 
wall reduces the build-up time of the blueshift, and thus 
would increase the final 1 -I- z. So, the difficulty here is in 
achieving a compromise between two contradictory goals: 
making the minimum of 1 -I- z on rays within the GRB 
sufficiently small, and the increase in 1 -|- z on rays within 
the afterglow sufficiently fast. 

As already stated in Sec. XHI, going over to the more 
general Szekeres metrics [28, 29], [7] may help in achiev¬ 
ing the collimation of the GRB flashes. This is not a 
particularly pressing problem (the collimation is only a 
hypothesis [6]), but employing the Szekeres model could 
also modify the results of Sec. XH in interesting ways. 
For example, an anisotropic BB hump might result in 
reduced angular size of the GRB generator when viewed 
from certain directions. 

Throughout this paper, an ambitious approach was 
taken, trying to explain the observed GRBs as if they 
all arose by blueshifting the relic radiation. This was 
to avoid the accusation that the author made his task 


process, their values given in (9.19) are much smaller than the 
first ones that produced 1 + z consistent with (9.5); they were 
{Ao,Bo, Bi) = (0.016, 0.01, 0.09). 
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unduly easy. Therefore, it was assumed that the lowest- 
frequency GRBs originate as the lowest-frequency emis¬ 
sion radiation of hydrogen atoms, and correspondingly 
for the highest frequencies. But it is possible that there 
are several mechanisms of producing the GRBs, and that 
blueshifting from the last scattering is only one of them. 
In that case, one can also consider shifting high emission 
frequencies to low GRB frequencies. For example, the 
blueshifts needed to shift the and the r'Hmax given 
by (9.3) and (9.2) to the minimum gamma-ray frequency 
given by (9.4) are, respectively, 1 -I- Zint ~ 9.1 x 10“"^ and 
l+2opt ~ 1-33 X 10“^. These values would be much easier 
to achieve than those considered in this paper. It would 
suffice to make the hump lower by reducing and i?o, 
and then suitably reducing Ate- The result would be as 
described in paragraph 2 of this section. 


XVII. CONCLUSIONS 

The cosmological model employed here is a spatially 
homogeneous Friedmann background (see Sec. VIII) into 
which a suitable number of Lemaitre - Tolman regions 
(see Secs. Ill and IV) is matched. Use was made of the 
long-known property of the L-T model that light rays 
emitted radially near those points of the Big Bang where 
the bang-time function is nonconstant (dts/dr ^ 0) dis¬ 
play arbitrarily large blueshifts instead of redshifts [3-5]. 
The present paper tested the hypothesis that the gamma- 
ray bursts (or at least some of them) arise by this mecha¬ 
nism. The blueshifts are generated only in the initial seg¬ 
ments of the rays, before they exit the extremum-redshift 
hypersurface (see Sec. VI). In the further part of the 
rays’ journey, redshifts are generated that can wipe out 
with excess the earlier-acquired blueshifts. The techni¬ 
cal problem to solve was this: Can the parameters of the 
model be chosen so that a substantial part of the blueshift 
survives the journey to the present observer and the ob¬ 
served frequency of the ray is in the gamma-ray range? 
Section IX answers this question in the positive. 

The models discussed here satisfactorily reproduce the 
GRB frequency range (property (I) in Sec. II, see Sec. 
IX) and account for the large distances to the GRB 
sources (property (5), see Sec. XIV). However, if the 
mechanism of generating the GRBs is such as discussed 
here, then the distances reported in the literature are 
calculated by an incorrect method and are heavily un¬ 
derestimated, see Sec. XIV. 

The other properties listed in Sec. II are accounted for 
only qualitatively, with varying degrees of success. The 
duration of the GRBs (Property (2)) comes out too long, 
but Model 3 of Sec. IX contains a free parameter that 
can control this quantity, see Sec. X A. To actually carry 
out the control, a much higher numerical accuracy would 
be necessary than was available to this author. 

The afterglows (Property (3)) necessarily exist in these 
models, but they last longer than observed and no way 
of controlling their duration was provided (see Sec. XI). 


As stated in the previous section, the observed durations 
are determined by the varying intensity of the radiation, 
which was not calculated in this paper. 

In Sec. XII, rays that are nonradial with respect to 
the L-T region were discussed. It was shown that the 
angular size of the GRB generator implied by Model 3 of 
Sec. IX is about twice as large as the resolution of the 
GRB detectors. So, with moderate improvements in the 
model that are anyway needed for other purposes (see 
Sec. XVI), the perturbations in the GMB radiation im¬ 
plied by our model will not be larger than those actually 
caused by the GRBs. 

The collimation of the GRBs into narrow jets (property 
(4)) appears to the observer as an illusion, as explained 
in Sec. XIII. This property is not directly implied by 
observations, it is just a hypothesis. However, it may be 
possible to deal with it when the quasi-spherical Szekeres 
(QSS) metric [28, 29], [7] is employed instead of L-T. 

In Sec. XVI the improvements needed in the model 
were discussed. 

Even if the models presented here prove to be unsat¬ 
isfactory explanations of the GRBs, they say new inter¬ 
esting things about the physics (in particular, optics) in 
the L-T models. 


Appendix A: The age of the Universe at last 
scattering vs. redshift 

There is a spurious discrepancy between the data on 
the age of the Universe at last scattering, tls, and the 
redshift of GMB. Namely, once k = 0 and Ha = 0.68 
have been specified for the AGDM model, the metric is 
uniquely defined, and then, given redshift, tls can be 
computed or vice versa. The value of the redshift of the 
GMB radiation is currently given as 2 ls = 1090 [18, 19], 
and then the calculation based on the AGDM metric (7.4) 
leads to tls = 4.77 x 10® y given by (7.16). On the other 
hand, the most-often cited value is tls = 3.8 x 10® y [33] 
(which, by the same method, would imply zls = 1269.3). 

The solution of this seeming contradiction is this: the 
AGDM model of (7.4) does not apply before the last scat¬ 
tering because it implies zero pressure via the Einstein 
equations. Eor that epoch, the pressure of radiation and 
matter cannot be neglected, and a more general model 
must be used, which, given zls = 1090 implies a tls value 
close to that given above [34]. The redshift at last scat¬ 
tering is calculated by sophisticated methods of particle 
physics - see Ref. [35] for a readable account. 


Appendix B: Blueshifting the black-body spectrum 

The maximum observed intensity of the GMB radia¬ 
tion by today is /max ~ 5 x 10“^® lU, where lU 
W/(cm^x sr x Hz), at the frequency J^max Ri 2 x 10^^ 
Hz [24]. Since T/zzmax is constant (Wien’s law), it fol- 
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lows from Planck’s formula I{iy) = x exp [ hi./(kr) ] -i 
that /max oc t'max^- The radiation frequency obeys 
i^emitted = i^received(1 + z), by the definition of z. The 
redshift between the instant of emission of the CMB radi¬ 
ation and the present time is 1 -I-zls ~ 1091 by (7.11), so 
the maximum intensity of the relic radiation at the time 
of emission must have been /max x (1-I-2ls)^ ~ 6.49 x 10“® 
lU. If this were blueshifted by 1 - 1-2 = 10^ to the gamma- 
ray frequency range, the resulting black-body radiation 
would have the temperature 3 x 10® K and peak inten¬ 
sity equal to 6.49 x 10® lU. This is many times more 
than observed ~ exemplary GRBs have intensities below 
10“®4 W/(cm® Hz) [36]. Thus, the GRBs cannot arise 
by blueshifting the thermal radiation with the black-body 
spectrum preserved “ and indeed the spectra of the GRBs 
are not thermal [21]. 


Appendix C: Remarks on choosing the shape of the 
BB hump 


the affine parameter A. With too large a step the inte¬ 
gration jumps over the LSH, and crosses the BB without 
noting that it had crossed the LSH in the same step. A 
very small step along the whole ray would make the cal¬ 
culation prohibitively long. Sometimes the same problem 
appeared in a vicinity of other locations, for example, in 
crossing the ERH or in determining the instant of obser¬ 
vation of a ray. In such cases, the step in A was changed 
in flight to either increase the precision near to critical 
locations, or relax it to accelerate the calculation. 

For example, in Model I, A A was 

AAi = 2.5 X 10"®, (Dl) 


all along the ray displayed in Fig. 2, which crossed the 
LSH in step # 31,184,648, with 


l + z 


LSHl 


9.11448916340873438X lO"!® 
-0.13941448281939275 NTU 


and it crossed the BB in step # 31,184,684 at 


(D2) 


As the first Ansatz for the shape of the BB hump, the 
Gauss-type family of curves was chosen 

tsir) = tBi + Ae"^’'^ (Cl) 

with A and B being adjustable constants and tef given by 
(8.3). However, the two free parameters did not provide 
sufficient flexibility. The rays (integrated backward in 
time) either hit the LSH with insufficient blueshift, or 
escaped back out through the ERH and acquired large 
redshift at the LSH. The strongest blueshifts were of the 
order of 1 -I- 2 « 10"® instead of the desired 10"® or less. 

The next shape that was tested was similar to that 
in Fig. 1, except that the upper left curve was also an 
ellipse. With this, the range of blueshifts given by (9.5) - 
(9.8) was achieved, but the humps were rather high and 
wide. In one model, which accounted for (9.7), the time 
difference between the maximum of ts and its flat part 
was 0.21T, with T given by (7.10), and the comoving 
radius of the hump was = 0.09, with the observer 
sitting at r ^ 0.39. In the other model the hump was 
only 3 times lower, and just as wide. 

The shapes of the humps that were finally used (see 
Sec. IX) emerged as the next corrections with respect 
to the one with two ellipses. They resulted in a radical 
reduction of both the height and width of the hump. 

Appendix D: Calculating 1 -|- 2 at the LSH 

In the neighbourhood where the rays from Figs. 2 and 
3 cross the LSH, it lies extremely close to the BB. In 
Model I the difference in t between them along the ray 
defined by (9.13) is ~ 5.54 x 10"^® NTU « 5.43 x 10"^ 
y ~ 20 days. Therefore, determining the instant when 
the ray crosses the LSH, and the value of 1 -|- 2 at that 
moment, requires a suitably small numerical step A A in 


tBBi = -0.13941448281944963 NTU (D3) 

with 1 -|- 2 at the level of 10"®®. The exact value was not 
captured by the program (the last 1 - 1-2 that was reported 
was negative, i.e. arose after crossing the BB). 

For Model 2, the increments in A were 

AA 2 = 2.5 X 10"®. (D4) 

The ray from Fig. 3 crossed the LSH in step # 5,175,194 
with 

/ 1 -k 2 \ _ / 1.07858890707746014 x 10"'^ \ 

it/ ~ \ -0.13945300368712291 NTU / ’ 

(D5) 

and it crossed the BB in step ff 5,175,204 with 

tBB 2 = -0.13945300368714261 NTU (D 6 ) 

and 1 -I- 2 ~ 5 X 10"®. The problem with capturing the 
last 1 - 1 - 2 , mentioned under (D3), recurred also here. 

Appendix E: Problems with calculating k'" from 
(4.13) 

Glose to the BB, where /? ~ 0, the expression under the 
square root in (4.13) (call it £ 1413 ) becomes a difference of 
two very large numbers. This follows from (4.10): where 
/? —>■ 0 with C ^ 0, necessarily (fc*)® —>■ 00 . The £1413 
must be positive, but sometimes the program could not 
ensure £1413 > 0 at small R. In calculating the paths 
of the rays in Figs. 10 and 12 - 14, whenever £1413 was 
negative, the program would replace it by [£1413 ], on the 
assumption that £1413 < 0 resulted from numerical errors 
in calculating £1413 ~ 0 . 

Along every ray in Figs. 12 - 14, fc’’ was expected to 
change sign from — to -I- at a certain point. In those 
cases, £1413 < 0 was a signal that fc’’ was near zero, and 
from that point on the program changed the sign of fc’'. 
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Appendix F: Deriving (12.6) 


Using this in (7.2) we find 


The equation of a radial null geodesic for (7.1), using 
(7.2) and r = 0 at the observer O, has the solution 

T] = rjo — In (^V—kr + \/l — . (FI) 

The value of rjo is calculated from (7.2), taking (8.3) for 
ts and t = 0 for the observer; it is 

r]o = 0.59249644326164175. (F2) 

With r = 0 assumed for the observer, r at D is the roa 
from (10.7), and then, from (FI) 

= 5.57539803067893525 x 10"^ (F3) 


to = ts + 


Mo 

(-fc)3/2 


(sinh?7D - rjo) 


= 1.14196789587062714 X 10"'^ NTU, 
Sd = 3.886639547232007 X 10-^ 


(F4) 

(F5) 


and from here (12.6) follows. 
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